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Abstract 

In this paper, thick iC-branes refer to four-dimensional domain walls generated by a scalar 
field with non-standard kinetic term. We first give a gauge independent description of the linear 
perturbations of a wide class of thick i^-brane models. Then we decompose the original metric 
perturbation into 4 scalars, 2 transverse vectors and a transverse-traceless (TT) tensor. We show 
that each type of these perturbation modes evolves independently. As compared to the standard 
case, equations of both vector and tensor perturbations are not modified. On the contrary, the non- 
standard kinetic terms usually cause modifications to the mass spectrum and effective potential of 
the scalar perturbation. In order to reduce the scalar perturbation equations into a Schrodinger 
like equation, a constraint on the matter Lagrangian is needed. This constraint might be useful in 
modeling ET-branes. 
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I. INTRODUCTION 



It is not a new idea that our world might be some topological defects, such as a Nielsen- 
Olesen vortex in a six- dimensional jlj], or a domain wall in a five- dimensional flat space- 
time. Instead of the conventional idea of "compactification" , a creative idea called dynamical 
localization was applied in these models. As a result, non-compacted extra dimensions 
are allowed, and in the low energy everything are trapped on these topological defects. 
In the 1990 's, it was found that warped extra dimensions can be very useful in solving 
some fundamental problems, such as the hierarchy problem and the cosmological constant 
problem 0, 0] , see also [3] for reviews. In 0, 0| , branes were called thin branes, because 
the thicknesses of branes were assumed to be zero. 

It is a natural idea to generalize the original vortex world or domain wall world scenarios 
to space-time with warped geometry [8-25]. Such models are also known as the thick brane 
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models, because branes in these models usually have non-zero widthes, see 
of thick brane solutions in various interesting cases. For some of the recent developments 
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One of the important issues in brane world models with large extra dimensions, is the 
trapping of various kinds of bulk matter fields on either thin or thick branes, see 39M45| for 
the localization of g auge fields, and 46N58| for fermion and scalar fields. The localization of 
j9-form fields 59h66|. and the Elko spinors 67j were considered recently. 

However, before we discuss the localization of bulk fields, we should first verify that 
the brane solution is stable. To discuss the stability of branes, we need to consider linear 
perturbations of the background brane solution. If no tachyon appears in the spectra of 
the perturbation modes, we say the background solution is stable, otherwise, the solution 
is unstable, see [l7J, |l9|. Besides, to reproduce the four- dimensional gravity on the brane, 
we need to localize the zero mode of tensor perturbation 0, 15, 17, 68]. Some modifications 
to four-dimensional Newtonian's law of g ravity can also be obtained by analyzing the linear 
structure of thick brane models 15|, |l7J, |69j- Thus, the linearization is very important in 
modeling thick brane models. 

In most thick brane models, branes are generated by single (or multi) real scalar field(s) 
with standard kinetic term(s). However, it is also possible to construct thick branes by using 



scalar fields with non-standard kinetic terms 70-81]. Scalar fields with non-standard kinetic 



terms (also called as K-fields) appears in the low energy effective action of string theory, and 
can be applied to describe the inflation of early Universe 82j-|84j. So the generalization from 
the standard thick branes to the iT-branes is both interesting and physically im por tant. 



The fermion localization on some K-branes, was considered, for example, in 72|, |78| . As 



to the issue of linearization of f^-brane, some previous works can be summarized as follows: 
stability of the domain wall solution under linear perturbation of the 
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scalar field was considered. However, since these models involve gravity, the metric 
perturbation should also be considered. Usually, the metric perturbation couples to 
the scalar field perturbation, which renders the problem rather difficult to solve. 



The tensor perturbation of thick K-branes and the localization of gravity were dis- 
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cussed in [76] . It was found that the tensor perturbation decouples from the K-G.e\d 
perturbation. 



The scalar metric perturbations of thick K-brane are, to our knowledge, not yet been 
discussed systematically. In fact, the issue of stability of thin brane models under scalar 
perturbations has been investigated in [85H87J by using non-standard bulk scalar fields. 



However the scalar fields in 85 



87] are not used to generate the brane, but to stabilize 



the distance between thin branes via the Goldberger-Wise mechanism 88j, |89(. So, 



the branes in 85M87|| are not i^-branes in our definition. Therefore, we need a refresh 
discussion on the scalar perturbations of thick iT-branes. 

In this paper, we try to give a general and systematical discussion on the linearization of 
a large class of iT-brane models. In section [Til we present the i^-brane model in a general 
form, and give the background dynamical equations. The linearization of our model as well 
as the issue of gauge freedoms in linear perturbation equations are generally discussed in 
section [TTT1 In section HV] we simplify the linear perturbation equations by using the scalar, 
tensor and vector (STV) decompositions of the metric perturbations. Some gauge-invariant 
variables are constructed and used in the perturbation equations. The issues of stability, 
localization of metric perturbation are discussed in section |V] In the end, we give a brief 
summary and stress some possible applications of our results. 



II. i^-BRANE MODEL AND BACKGROUND EQUATIONS 

In the model with a single K-field, the most general action for thick branes is given by 

S = S 9 + S m = J d 5 xy/=g~ f-^R + £(0, X)j , (1) 

where k\ = 871G5 is the gravitational coupling constant and G5 is the five-dimensional 
Newtonian constant. We use X = — \g MN V n4> to denote the kinetic term of the 
scalar field, so that the standard lagrangian density is given by C = X — V(<p). In this 
paper, we always use fi, v = 0,1,2,3 to denote the indices of brane coordinates, and use 
M, iV = 0,l,2,3,5to represent the indices of bulk coordinates. 
The Einstein equations are 

Rmn — ~9mnR = I^TmNi (2) 



where 



T MN = — = j-^ = g MN C + £xVm0Vjv0 (3) 



is the energy-momentum tensor. 

The metric for domain wall brane is assumed to be 



ds 2 = a 2 (y)r] tMJ dx^dx" + dy 2 , (4) 
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where rj^ = diag(— , + ,+ ,+) is the four- dimensional Minkowski metric, y 
the extra dimension. The warp factor a 2 (y) is an even function of y. 
It is more convenient for us to start with a conformal flat metric: 



= x 5 represents 



9mn = a 2 (r)rj MN . (5) 

This metric relates with metric fll|) via a redefinition of the fifth coordinate dr = a~ 1 dy. 
Starting from the conformal metric, we immediately obtain the expressions for the connec- 
tion, the Ricci curvature and the scalar curvature 

^m N = \s^d M a + 5^d N a - r] MN r] PQ d Q a), (6) 

D d M ad N a d M d N a ( a'\ 2 a" 

Rmn = ® 5 3 2rjMN — -Vmn~, {') 

a z a \ a J a 

correspondingly. Primes denote derivatives with respect to r. Plugging these expressions 
into eq. (j2J) and use eq. (jBJ) we obtain the Einstein equations 

d M ad N a n d M d N a a" 

6 3 V Srj MN — 

a a a 

= i4a 2 r) MN C + K 2 5 C x y WV ' N (j), (9) 

or in component form: 

6^-^ = K 2 5 a 2 C + K 2 5 C x <f>' 2 , (10a) 



3^ = 4a 2 C. (10b) 
The equation of motion (EoM) for the scalar field is 

£^ + V A/ (L x V M 0) = 0. (11) 
As our metric (JSj) is considered, we have 

<f>'C' x + C x U" + 3^/) = -a 2 ^. (12) 

This equation can be derived from the Einstein equations as a natural result of the Bianchi 
identity V n Gmn = 0. Thus, the Einstein equations ffTUj) constitute the independent and 
complete background dynamical equations and will be used frequently in the following dis- 
cussions. 
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III. LINEARIZATION AND GAUGE FREEDOMS 



In principle, one can find a lot of thick brane solutions by solving equations f llUp and 
(|12p . However, not all the solutions of the dynamical equations are stable. To discuss 
the stability of thick .fT-branes, we have to analyze the linear structure of our model. If 
some tachyon modes appear in the spectrum of linear perturbations, then the corresponding 
i^-brane solution is not stable. 

In our model, K-fie\d couples with gravity, so perturbations from both the metric and 
the i^-field should be considered, so 

9mn = 9mn + Sg MN , = + 50, (13) 

where Qmn and are the background solutions which satisfy equations ffTUl) and f|T2|) . SgMN 
and 50 are small perturbations around the background solutions, both of them are functions 
of the whole bulk coordinates x M . It is convenient for us to redefine the metric perturbation 

as 

SgMN = a 2 h MN (x p ,r). (14) 

Since the perturbed metric gMN should also satisfies the orthogonal relation gupg PN = S M , 
we immediately obtain 

Sg MN = -a 2 V MP V NQ h PQ . (15) 
Before we keep on the next step, let us first state some notations: 

1. Note that usually d p = g p ®dQ. However, since we have used the conformal flat 
metric, we frequently encounter term like f] P( ^dQ. So it is more convenient for us to 
denote d p = r] p< ^dQ and d p = rf v d v for simplicity. For the same reason, we denote 

h MN = rjMPjjNQhp^ so that 5g MN = _ a 2 h MN and h = jjMN^^ 

2. = r] MN d M d N = d M d M and CK 4) = rj^d^ = are d'Alamber operators 
defined with Minkowski metrics and the ordinary partial derivatives, obviously, D*^ = 

+ d r d r . 

3. The symmetrization of the indices M, N of a second order tensor Tmn is realized by 
using the symmetrization bracket, which is defied as 

T(mn) = 77 {Tmn + Tnm) ■ (16) 



Note that the perturbation of the Einstein equation ([2]) is 

1 „ „ 1 



SRmn — 2^9mnR — ^9mnSR — k 5 STmn- (17) 
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Plugging eq. flS}, and fjA2[) - fjA5H into above equation and after a simplification, we obtain 

d P d (M h N)P - \ ( □<» + 3 a -d r ) h MN - 2 ± d J^> 
2 \ a J a a 



a' fa'\ 2 1 

+3—d( M h N ) r + 21 — 1 r] r{M h N ) r - -d M d N h + r] MN X 
Ob \ a j £ 



-CxxaT 2 (f>'S<t>' + ^£xxa-y 2 /w) , (18) 



where 



X = l n W h _ id p d Q hpn + --«9 r /i - 3-d p /i Pr 
2 2 g 2 a a 

a" , , , 1 



- 3— /i rr - ^(^£ x /2 /i rr - £ X <W + a 2 C^4>). (19) 

In component form, the perturbation equations are 

1 3 a' 

(//, i/) : d^d p h u)P - -D^h/ju, - -—fyk^ 

1 a' 

- -d^dyh + 3—d(nh v)r + r^X = 0, (20a) 

2 a 

(//, r) : d^d p h r)P - in( 5 ) V + Ifahrr - \dAh = i^Cxtfd^, (20b) 

1 3 ft' 1 

(r, r) : d r d F h Pr — -D^h rr + ~—d r h rr — -d r d r h + X = 2nlCx4>'S(p' 

+4<P' 2 (C x ^(f) - C X xa- 2 <j)'S<J)' + \cxxa- 7 ^ 2 h rr ). (20c) 

These equations describe the evolution of the linear perturbations for the most general case. 
However, there are some gauge freedoms still exist in the perturbations. To see why, let us 
recall the general covariant principle, which states that the dynamical equations should be 
independent from our preference of coordinate systems. So under an infinitesimal coordinate 
transformation 

x M ^x M = x M +£ M (x p ), (21) 
an arbitrary second order tensor should transforms as 

Tmv(^) = §M^r JK (x% (22) 

Just like we have done in eq. ( JT3]) . we can split Tm n as a background field plus a small 
perturbation: 

Tmn(x q ) = Tmn(x q ) + 5f MN (x Q ), (23) 
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then, under linear approximation, transformation in ([22]) is equivalent to the following gauge 
transformation of the linear perturbation: 

6Tmn(x q ) = 6Tmn(x q ) - d M t P T PN (x Q ) 

- d N i P T PM (x Q ) - ed P T M N{x Q ). (24) 

As the metric tensor is concerned, we have 



Ah M N = h 



MN 



(25) 



where £m = Vmn£, n and "A" indicates the change of perturbations. 

For scalar field perturbation, the gauge transformation can be obtained by following 
similar discussions (see |90j for details): 

A5<p = -C ft. (26) 

It can be shown that, under the gauge transformation (125]) and ( 126]) . the linear perturbation 
equation (ITS]) is invariant. In other words, (5<p, Iimn) and (8(f), Iimn) satisfy same equations. 
This is an important information, because we can use these gauge freedoms to eliminate 
some of the perturbations. A simple example is that we will have 5<p = 0, if we take 



(27) 



Likewise, we can eliminate some other perturbations by using the residual freedoms of choos- 
ing 

So, on one hand, we have to eliminate all these freedoms because they are not physical. 
While on the other hand, these gauge freedoms allows us to working in some particular 
coordinate systems in which the perturbation equations are relatively easier to solve. Indeed, 
we have different ways in eliminating the gauge freedoms. 

The first choice is to take gauges directly, just as we have shown in the example above, 



see also [9ll-l94|. Some authors prefer to take gauges in the light-cone coordinates |95N99 



However, it is difficult to eliminate all the gauge freedoms completely if we directly take 
gauges. 

Nevertheless, if the space-time is highly symmetric, we can take gauges after decompose 
the metric perturbation into scalar, vecto r and ten sor parts (the STV decompositi on) . This 
method was first used in cosmology, see 100Ml03| ] for original contributions and 104H107 



for recent applications of this method in cosmology. The advantage of this method is that 
we can not only find some gauges which completely eliminating all the gauge freedoms, but 
also can construct gauge invariant quantities which can be treated as the physical dynamic 
variables when we try to quan tize the p erturbations 103]. This method was also generalized 
in brane world models 17l. ll9l . ll08l - lll5j to analyze the stability of thin or thick brane models. 



Beside these methods, the covariant and gauge- invariant approach is also a widely used 



method in analyzing perturbations 116| . This method is based on the 3 + 1 deco mpos ition 
of the space-time. For applications of this approach in brane world cosmology, see [ll71 | and 
references therein. 

In this paper, we wi ll fo llow the second approach, the details can be found in standard 



text books such as |90|, Il06 



7 



IV. STV DECOMPOSITION AND GAUGE INVARIANT VARIABLES 



It is always possible to make the following decomposition (see 106| for how to realize 
this idea in cosmology): 

V = dfj,F + Gfj,, (28a) 

V = %vA + d&B + 2d {fM C u) + D^, (28b) 
where C^, G M are transverse vector perturbations: 

d»C, = = d"G„, (29) 
and Dfj, v is transverse and traceless (TT) perturbation: 

VD^ = = D». (30) 

Here all indices are raised with rf" . 

Using the properties of the decomposed metric perturbations, we can now figure out how 
these perturbations transform under gauge transformation (|25|) . For scalar perturbations, 
we have 

AA = -2-C, Ah rr = -2C - 2-C, (31) 
a a 

AB = -2C, AF = -C - C'. (32) 

We have used the decomposition £ M = + with <9 M £ ±M = 0. 
For vector and tensor perturbations we obtain 

AC7 M =-#, AG, = -^', (33) 

and 

AD^ W = 0, (34) 

respectively. It is useful to define ip = F — ^B' and v^ = G^ — CJ_ . Obviously, = — £ r , 
while Vp, D^ u and the following scalar quantities are invariant under gauge transformations: 

S = h rr — 2— (aip)' , 
a 

= A-2-ifj, 
a 

$ = 6<p - (p'ip. (35) 

Now let us discuss the perturbation equations. Plugging eqs. f l28|) into eqs. fl20|) . and 
using expressions in (1A6j) we obtain the following equations: 

a' 1 

if)' + 3—ib - A - -h rr = 0, (36) 
a 2 



X=-U^A+-A" + --A' 
2 2 

V") + d {lx v' v) 



_ I', (37) 

2 2 2a y 1 

3— d^v v ) + dfrv' v) = 0, (38) 

+ 2^ + 3^ = 0, (39) 
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from the (//, v) components. 

By using ffT9]) and (IA6D . one can immediately obtain 

1 = in( 4 M + 2A" - --hrr' + 6— A 1 - 3—h rr 

2 2 a a a 

- ^Cx^Kr - C X <P'H' + a 2 L^0), (40) 

However, as we will show below, we need not to plug this complex expression into the 
perturbation equations (for example, eq. (|37|) ). 
From the (/i, r) components we obtain 

□ (4) *V = 0, (41) 
—Kr-^A^^Cxm- (42) 
Finally, the (r, r) component can be written as 

D (4y _ \u^h Tr + --h„' - 2A" + X = 2k\L x 4>'^ 

W^Cx+St - C xx ar 2 <l>'6<l>' + l -C X xa- 2 (f>' 2 h rr ). (43) 

Obviously, different types of perturbation modes evolve independently. So we can analyze 
them separately. 

The equations for vector and tensor perturbations are obviously gauge invariant. One 
can check that all the equations for scalar perturbations are also gauge invariant. If we plug 
the gauge invariant quantities defined in eq. (155]) into the scalar perturbation equations ( 136]) 
and ( 142]) . we would get 

-¥ - \E = 0, (44) 
^Z-h> = 4Cx<t>% (45) 
The third independent equation for scalar perturbation comes from eq. ( 143]) 

_ h" - - ^CxxaT 2 ^ 

2 2 2 a 5r ^ 

= 2k 2 £ x 0'$' - Kl^Cxxory^' + k\^L x ^. (46) 

Here, to eliminate X in eq. ( 143]) . we used eq. ( 137]) rather than eq. ( 140]) . 

Using gauge-invariant variables is one of the natural ways to eliminate the gauge freedoms 
in the perturbation equations, another way to achieve this goal is to take gauges. For 
example, by taking tp — 0, B — and = 0, we can eliminate the gauge freedoms in £ r , ( 
and ^ J " At , correspondingly. Under this gauge, {H, $, v^} simply reduce to {h rr , A, (/), G M } 
respectively. So, the final equations for {h rr , A, <fi} are nothing but eqs. (144]) -( l46]) . The 
gauge we used here is called the longitudinal gauge, or conformal Newtonian gauge, which 
is widely used in literatures. Therefore, taking longitude gauge is totally equivalent to using 
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the gauge invariant variables defined in eq. (135]) . Both methods completely eliminate the 
gauge freedoms in the perturbation equations. 

Note that the construction of gauge-invariant variables is not unique. In fact, for any 
gauge that completely eliminates the gauge freedoms, we can always construct the corre- 
sponding gauge-invariant variables. A good choice of gauge usually helps us to simplify 
the perturbation equations. Usually, we can simplify the perturbation equations by simply 
taking the longitude gauge. However, for some special models, such as thick branes in f(R) 



gravity 32], there might be some other gauges which are superior than the longitude gauge. 



V. STABILITY AND LOCALIZATION OF ZERO MODES 

Let us first consider the vector perturbation which is described by eq. (138]) and eq. (|4TI) . 
If we expand as 

v^e.e^xir), 5% = 0, (47) 

then, eq. f|4Tl) indicates that m 2 = —p\p x has to be vanished, namely, only the zero vec- 
tor mode survives. Thus, any K-brane solution in our model is stable under the vector 
perturbation. From eq. (1381 . we get the expression of the vector zero mode: 

Xo(r) oc a' 3 . (48) 
We can also expand tensor perturbation described by eq. (139|) . as follows: 

D, v = V«" 3/2 e ip ^ A p(r), d»d%„ = = /TV (49) 
Then, eq. ( 13T?]) can be written as (see also [761 ]) 

d r + 1-) (d r - 1-) p(r) = -m 2 p(r). (50) 



2 a J V 2 a, 

In the theory of supersymmetric quantum mechanics, the above factorized equation indicates 
that m 2 > 0, and thus fT-brane solutions satisfying eqs. flTUj) are also stable under tensor 
perturbation. Note that the tensor zero mode is 

p (r)oca 3 / 2 . (51) 

In order to localize the zero mode of vector or tensor perturbation, we need 

dr\ Xo (r)\ 2 = l or, / dr\p (r)\ 2 = 1. (52) 



However, as stated in [17l |. these two requirements cannot be satisfied simultaneously. In 
addition, if we demand a finite four-dimensional Planck mass, the zero mode of tensor 
perturbation should be localized on the brane. As a result, the vector zero mode cannot be 
localized. 

Since the introducing of K-field does not change the equations for both vector and tensor 
perturbations. There is nothing new in the results as compared to the standard case C = 
X — V(<fi). Now, let us consider the non-trivial case, the scalar perturbations. 
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Plugging eqs. (144]) and (145]) into eq. (146]) and after simplification, we obtain 



(53) 



where 



01 = # ln(£ x 0') - 4^ - 2^0' 
2^-H-d r \n(C x ^) 



-X- 



X 



XX 



£ 



X 



a 2 = 2d r \n(£ x <p') - 3H - 



£ 



x 



-X 



X 



xx 



«3 



2 C x 



-d r \n(Cx(p') + H 



(54) 

(55) 
(56) 



In order to analyze the mass spectrum of scalar perturbation, we hope to transfer the eq. ( 153]) 
into a Schrodinger like equation. So we hope «3 is a constant. The general form of £ that 
satisfies this constraint is 

£ = 1 f\4>)X n - V(cP). (57) 
Here /(</>) and V(</>) are arbitrary functions of 0; 7 and n are real numbers, so that 



£ 



xx 



£ 



X 



and, 



X = n-1, 



J 1 



£ 



£ 



X 



f 



(58) 



ax = 2H^- - 2(n + l)H' + (n - 1)# 2 , 



a 2 = d r In 



J 



7 n+2 



«3 



(59) 



Here we have defined J = (£x<f)')^ r ^ f~ x and H = — . This Lagrange density contains a lot 
of interesting models. 
Expanding \&(r, x M ) as 



ty(r, x A 



7 n+2 



we can transfer eq. ( 1531) into the following Schrodinger equation: 

m 2 y?(r) = -ip(r)" + V cS (r)(p(r), 

where 

2m 2 



m 2 



l + n)' 



(60) 



(61) 



(62) 
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and the effective potential reads 
V eS (r] 

the coefficients are given by 

ci 

c 3 



a' J' 



ci — + c 2 7 + c 3 — 



a J 



J' 



J 



a" 1 J" 
a 2 J 



(63) 



(4 + 9n + 6n 2 ) 

(1 + n) 2 
(2 + n) 



2n 



c 2 



1 + n 



2 ' 



1 + n* 



2 ' 



c 4 



(2 + 3n) 
(1 + n) ' 



(64) 



For the canonical scalar field, C = X — V, we have f((f>) : 
The coefficients of the potential is c\ = ^, c 2 = f , c 3 



Kff(r) 



2 a A\a) a (j)' ) cf)' 



1, n 


= 1 and J = 


3 

4' 


c 4 = |, and 


(8 






0" 



y 2 , m 2 



m 



(65) 



which consists with the results in 17, 19, 108, 110 



For general case, both the mass spectrum and the effective potential are modified. When 
n < — 1, we get m 2 < 0, so any solution is unstable in this case. 

Obviously, the effective potential is determined only by the background solution. Thus, 
as a particular i^-brane solution is given, we can discuss its stability by simply analyzing 
the corresponding effective potential. But we leave such special discussions to our future 
works. 



VI. CONCLUSIONS AND OUTLOOKS 

In this paper, we generally discussed the linearization of at thick X-brane models. We 
obtained the most general linear perturbation equations which are independent from any 
gauge choice. Then we applied the scalar, tensor and vector decomposition to the metric 
perturbation. It turns out that each type of the perturbations evolves separately. Then, 
by defining a group of gauge-invariant variables, we rewritten the perturbation equations 
into some Schrodinger like equations. Note that in order to write the scalar perturbation 
equations into a Schrodinger like equation, we imposed a constraint on the Lagrangian of 
scalar field, namely, 

£ = 1 f(<P)X n -V(<j ) ), (66) 

where f(<f) and V{4>) are arbitrary functions of <f>; 7 and n are real numbers. For thick 
brane models with such a matter Lagrangian, the scalar perturbation equations can always 
be recast into a Schrodinger equation. 

There is no modification for tensor and vector perturbation spectra in thick X-brane 
models, no matter what the matter Lagrangian is. Thus, as to the standard case C = 
X — V((f>), if we want to have a finite four- dimensional Planck mass, we have to localize the 
tensor zero mode. Then, the vector zero mode cannot be localized. On the contrary, the 
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appearance of non-standard kinetic terms does impact the mass spectrum and the effective 
potential of the scalar perturbation. We will discuss applications of our results of scalar 
perturbations to particular models in our future works. 

We would like to stress that by using the longitude gauge, we cannot write the scalar 
perturbation equations of some models, such as domain wall branes generated by tachyonic 



fields |70|, [79|, into a Schrodinger equation. But, this problem might be solved by working 
in some other gauge choices. Since we have obtained the gauge invariant perturbation equa- 
tions (|36|) - (|43|) . we are ready to take any appropriate gauge or construct the corresponding 
gauge-invariant variables. We will consider this possibility in our following works. 
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Appendix A: Appendix 

The linear perturbations of the connection, Ricci tensor, Ricci scalar are given by 

S^mn = 9(Mh^) - h MN -d p h MN + r)MN—hr, (Al) 

; a 2 a 

5Rmn = dpd^jh^ — -D' 5 ^MAf — 2~a^ r ^ MN 

huN ~ 2 ( — ) h>MN + VMN—dph^ 

a \ a J a 

Jl / i\ 2 



+VMN — Kr + %VMN ( — J Kr - ^ 9 M 8 N h 



a \ a J 2 

1 a 'p)h^^ a 'f) h q q/ d (Mah N)r 
— -i]MN—o r ri + o—Of M n N ) r — 2 

2 a a a a 

+2 r] r{M h N ) r , (A2) 



a 2 SR = ( d M d N h MN - n®h - i-d r h + 8-dph? 
\ a a 



a" . I a 



+8—h rr + 4 ( - ) h rr ) , (A3) 



respectively. 
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The perturbation for the kinetic term X and the energy-momentum tensor Tmn are given 
as follows: 

5X = a- 2 (hi rr (j)' 2 - <j/6<f/) } (A4) 
ST MN = a 2 h MN C + V m0V ' N (j)(L x <t>&4> - C X xa~ 2 <P '&4> 
+^£xxa- 2 (f)' 2 h r ^ + 2C X V {M 4>V N) 5<P 

+Vmn (^Cx4>' 2 hrr ~ Cx<j>'S<f>' + a 2 £^) . (A5) 

Using the STV decomposition (|28|) . we get 

h = h rr + AA + n {4) B, (A6a) 
d M h Mr = d r h„ + D (4) F, (A6b) 
d M d N h MN = d r d r h rr + n<$(A + D^B + 2F'). (A6c) 
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